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ABSTRACT 


In the Haystack Pointing system, errors caused by atmospheric refraction, 
gravitational deformation of the antenna, skewed axes and resolver error are com- 
pensated for by a correction program in the computer which adds the necessary 
biases to the geometric values of azimuth and elevation to produce an effectively 
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HAYSTACK POINTING SYSTEM: RADAR COORDINATE CORRECTION 


be INTRODUCTION 


The Radar Coordinate Correction program of the Haystack Pointing system 
accepts as input a pair of angles (8, $) which are the azimuth and elevation coordin- 
ates of a point at which the antenna is to be pointed. The program modifies these 
angles by adding corrections derived from tables using the method of table look-up and 
interpolation and produces as output a pair of antenna pointing angles (6. o.) which 
are corrected for gravitational deformation of the antenna, atmospheric refraction, 


resolver error, and skew of axes. 
II. PROGRAM OPERATION* 


The program is composed of two parts, the initializing package and the worker 
package. Entrance to each package from a calling routine is made with an RJP 
instruction. 

The initializing package sets up standard correction tables during system 
initialization. If the correction program is called via the attention symbol route for 
reinitialization, the initializing package asks the operator to enter values for the 
various parameters used in computing the correction tables. See Fig. | for typical 
sequence of questions and answers. A vacuous answer (the carriage return alone) 
results in use of a prestored "standard" value for that particular parameter. The 
operator also has the option of omitting from use any or all correction tables. After 
all parameters have been entered, the package computes the correction values and 
stores them in the various tables. Control is returned to the calling routine from the 
initializing package with an EXIT instruction. 

If temperature or pressure should vary significantly during an experiment, 
these values can be re-entered and a new refraction correction table computed with - 
Out interrupting the experiment. ‘The time necessary to recompute the refraction 
table after the parameters have been entered via the console is approximately 4 
milliseconds. 

There are three tables involved in the correction of the pointing angles (see 
Appendix B). The table REFRACTBL corrects for atmospheric refraction and is a 


*Use of this document assumes knowledge of TN-1966-10, ''Haystack Pointing System: 
Control Structure, " by J. D. Drinan and A. A. Mathiasen. 


CORRECTION PGM 
INCLUDE REFRACTION TABLE (Y OR N) 
Y* 


USE STANDARD VALUES (Y OR N) 
N* 


ENTER TEMP IN DEG C 
p> FD, Ges 


ENTER ATM PRESSURE IN MB 
», ©.6, G8 a 


ENTER PARTIAL PRESSURE IN MB 
XX. X* 


INCLUDE AZ TABLE (Y OR N) 
y* 


INCLUDE EL TABLE (Y OR N) 
N* 


Fig. l. Typical question-answer sequence. 


function of elevation only. The argument varies from 0° to 90° and is more 
densely packed for angles below 25°. The parameters associated with this table 
are temperature in degrees centigrade and atmospheric (i.e. , total) and partial 
pressures in millibars. 

The equations used in computing atmospheric refraction corrections are based 
on the results of a paper published by W. R. Iliff and J. M. Holt. ! They verified by 
experiment that the total atmospheric refraction correction angle fr could be 


accurately estimated by an equation of the form 


where N. is the surface refractivity, a g? b é are functions of the elevation angle 
rw. R. iff, and J. M. Holt, "Use of Surface Refractivity in the Empirical 
Prediction of Total Atmospheric Refraction," J. Research NBS 67D (Radio Prop. ) 
No. 1 (January-February, 1963). 


@. The value of the surface refractivity is computed using the Smith and Weintraub 


—. 
equation 


77.6 4810e 
© Peel Y Prarie |: 


where e is the partial pressure of water vapor in millibars, p is the total pressure 
in millibars, and T is the temperature in degrees centigrade. The values of a $ 
and ae will be precomputed and stored as constants in two tables called ATBL and 
BTBL for values of ¢ ranging from 0° to 90°. See Appendix A for calculation of 
ile and ae. 

The tables AZTBL and ELTBL are each functions of azimuth and elevation. 
They correct for antenna sag, resolver error, and skew of axes, with AZTBL con- 
taining the corrections for azimuth and ELTBL the corrections for elevation. Each 
table has associated with it a pair of vectors containing the arguments of azimuth 
and elevation. They range from 0° to 90° in elevation and 0° to 360° in azimuth. 
The arguments can be irregularly spaced. See Appendix B for the table format. 

The worker package uses the angle @ ina table look-up in REFRACTBL to 
produce a correction rT for elevation due to refraction. This correction is added 
to @ to produce the angle ¢@. The angles (9, ¢) are used as arguments in separ- 
ate table look-ups of AZTBL and ELTBL to produce correction values 65 and 53 
for azimuth and elevation, respectively. These values are added to (9, ¢) to pro- 
duce final corrected angles (6. $.): While the values of (9, ¢) are in revolutions 
(with a binary point at B27) the table look-up and interpolation are done in degrees 
with the correction values converted to revolutions (B27) just prior to being added to 
(6, ¢). Control is returned to the calling routine from the correction package with an 
EXIT instruction. The time required for one pass through the worker package is of 


the order of 2 milliseconds. 
Iii. TABLE LOOK UP AND INTERPOLATION 


The correction value for each table is computed according to the following 


equations: 


2. E. K. Smith and S. Weintraub, "The Constants in the Equation for Atmospheric 
Refractive Index at Radio Frequencies, "Proc. IRE 41, No. 8, 1035-1037 (1953). 


A. 8, 


The values stored in the refraction table REFRACTBLE are T.. 


The values stored in the elevation argument table REFRACARG are S- 


) — _ 
6, = 7+ Bi oo 5 Tit Th 
where sis o< Sia] 
B Oo 


The values stored in the azimuth correction table AZTBL are 6. ; 


The values stored in the azimuth argument table AZTBLAZARG are a 


The values stored in the elevation argument table AZTBLELARG are Ve 


or 

D..... "= §. .+ = {6. .,., —&) 

13] i,j ae i, tl i 
°= ) 

6s = De . + ; = ‘ 

2 i,j Vii-Y; aks | i,j 

where ene 
y= Oo aed 
C 5. 


The values stored in the elevation correction table ELTBL are € : 
The values stored in the azimuth argument table ELTBLAZARG are xy 


The values stored in the elevation argument table ELTBLELARG are 


2. * 62 ¢e——46 3) 
ie ist) a4 Rie 1 
6 er 
3 inj Yin) -Yj itl,j i,p , 
where x, <6@< Kael : 
0 =? “a 


For an elevation angle @ greater than 90°, the correction value by is 
computed using the angle 180° — @. The correction value 5. so obtained is sub- 
tracted from the angle ¢. Fora negative angle #¢, 5, is computed by means of 
extrapolation and is added to the elevation angle @¢. 

Interpolation is linear in both elevation and azimuth with an accuracy com- 
patible with the accuracy of the input angles. All operations are carried out in single 


precision, fixed point arithmetic. 


APPENDIX A 


The following equations were used in precomputing the values of ye and yr; 
which are stored as constants in the correction program. All equations and parameter 
values were taken from the paper cited in Reference lL. 


The refraction correction angle 7, measured in revolutions, is given by 


¢o Ss @ ~ 
where No is the surface refractivity. The values of a and De are given by 
a, = =e 
? (a + B) 


where @ is the elevation angle measured in degrees. The functions ay and a 


have the dimensions of degrees and degrees per surface refractivity unit, respect- 


ively. As Q@ becomes large, 7 approaches zero and b, approaches the product 


p 
of a constant times the cotangent of the elevation angle. The parameters A, B, C, 
D, E, F are positive constants and were determined empirically. 


Their values are as follows: 


A = 40.0 D = 42.5 
B = 2.7 E= 0.4 
C= 4.0 F = 2.64 


APPENDIX B 


The following discusses the contents and formats of the tables used in the cor- 
rection program. In all tables, the subscripts (i,j) vary with elevation and azimuth, 


respectively. 


REFRACARG REFRACTBL 


The elevation argument S. is stored as degrees with the binary point at B20. 
The values of a; and b. will be precomputed and stored in these tables as constants, 
with a binary point at B20. They are functions of elevation angle and are used in com- 
puting the refraction correction table REFRACTBL. The refraction correction value 
T. is computed as a function of temperature, atmospheric pressure, vapor pressure, 


i 
a. and b.. The value T is expressed in degrees with the binary point at B20. 


AZTBLAZARG 


AZTBLELARG 


ELTABLAZARG 


ELTBLELARG 


The arguments u., Vis Xi, y, are stored as degrees with the binary point at 


B20. The correction values 2 and G will be precomputed and stored in these 


> 


tables as constants, where 6. is the correction value for azimuth and G; . is the 
correction value for elevation. They represent the corrections necessary to account 
for errors other than the refraction error. Both 6. and | 5 are expressed in 


degrees with the binary point at B20. 
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